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Abstract 

The microscopic theory of current carrying states in the balhstic 
superconducting microchannel is presented. The effects of the contact 
length L on the Josephson current are investigated. For the temper- 
atures T close to the critical temperature Tc the problem is treated 
selfconsistently, with taking into account the distribution of the order 
parameter A (r) inside the contact. The closed integral equation for 
A in strongly inhomogeneous microcontact geometry {L^^q, is the 
coherence length at T = ) replaces the differential Ginzburg-Landau 
equation. The critical current IdL) is expressed in terms of solution 
of this integral equation. The limiting cases of L <C ^0 and L ^ S^q are 
considered. With increasing length L the critical current decreases, 
although the ballistic Sharvin resistance of the contact remains the 
same as at L = 0. For ultra short channels with L^od (a/j vf/(^d 
, is the Debye frequency) the corrections to the value of critical 
current Ic{L = 0) are sensitive to the strong coupling effects. 



1 



1 Introduction 



Weak superconducting links |jT| include the tunnel structures S — I — S 
(superconductor-insulator-superconductor) and the contacts with direct con- 
ductivity, S — N — S (iV is the normal layer) and S — c — S (c is a geometrical 
constriction) . Superconducting constrictions can be modelled as the orifice 
with diameter d in a inpenetrable sheet for electrons between two supercon- 
ducting half spaces (point contact), or as a narrow channel with length L in 
a contact with superconducting banks (microbridge) . Aslamazov and Larkin 
[0] have shown on the basis of a solution of the Ginzburg-Landau (GL) equa- 
tions that in the dirty limit and for small constriction's sizes L, (i -C ^(T) ( 
^(T) is GL coherence length) the S — C — S contact can be described by a 
Josephson model with the current-phase relation 

I = hsmip,h = nAl{T)/iAeRNn), (1) 

where Ic is the Josephson critical current, Aq the absolute value of the order 
parameter in the bulk banks, the critical temperature and the normal- 
state resistance of dirty microbridge. The critical current of microbridge (|^) 
depends on the bridge length as Jc ~ l/L. The expression (|I]) is valid within 
the domain of applicability of GL approach, i.e. for temperatures T close 
to Tc and L,d ^ C,o (^o — vp/Tc is the coherence length at T = 0, vp is the 
Fermi velocity). 

The present level of technology have made it possible to study the ultra- 
small Josephson weak links with the dimensions up to interatomic size. For 
example, it can be nanosized microchannels produced by means of a scanning 
tunneling microscope ^ or point contacts and microchannels obtained by us- 
ing the mechanical controllable break technique 0,1^,^. The microchannels 
between two superconductors can also arise spontaneously as microshorts in 
tunnel junctions[0], with the length L determined by the thickness of an in- 
sulator layer. The value of the critical current Ic of such microshorts presents 
the special interest in the case of tunnel structures based on high-Tc metal- 
loxide compounds. The case of small microconstrictions with dimensions of 
the order or smaller than coherence length , when the expression (|I|) for 
the critical current Jc ~ l/L is not valid, requires the microscopic consider- 
ation even for T near T^. Such microscopic theory of stationary Josephson 
effect in microconstrictions was developed in Ref. p[ for the ballistic channel 
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of zero length L = in the model of the orifice of diameter d << ^q. The 
Josephson current in this case is given by: 



7rAo(r) 
eRo 



sin if/ 2 tanh 



Ao(r) cos(/p/2 
2T 



TT < if < n, 



(2) 



(3) 



where S = nd'^/A is the contact cross-section area, iV(0) = mpp/i^Ti'^) is 
electron density of states at the Fermi surface. At temperatures T^ — T-^Tc 
expression coincides with the Aslamazov-Larkin result, Eq. (|l]), in which 
instead of the normal resistance Rn for dirty metal, the ballistic Sharvin 
resistance Rq (j^) is substituted. 

The purpose of this work is to present a microscopic theory of current 
carrying states in the ballistic microbridges having the length L arbitrary in 
the scale of the coherence length ^q- We have investigated the dependence 
of Josephson critical current on the ratio L/^q and analyzed the transition 
from the case of Ic{L = 0) to Jc ~ l/L (|l]) with increasing of the length 
L. 

In Sec. 2 we formulate the model of a microbridge and the microscopic 
equations for Green's functions with boundary conditions at the bridge edges. 
When the effects on the critical current of the microconstriction's length are 
studied, the crucial point, as always in space inhomogeneous superconducting 
state, consists in the self-consistent treatment of the order parameter distri- 
bution A(r) inside the weak link. In Sec. 3 the closed integral equation for 
the order parameter A in microchannel is derived for temperatures near T^, 
which in strongly inhomogeneous (L ~ ^o) microcontact geometry replaces 
the differential GL equation. The critical current Ic{L) is expressed in terms 
of the solution of this integral equation. The limiting cases oi L <^ C,o and 
^ ^0 are considered. We shall show that beside the characteristic scale 
^0 the length ajj ~ vp/ujd {ud is Debye frequency) appears in the case of 
ultra small channel. The length L ~ a/) is the length at which the frequency 
of the ballistic flight of electron from one bank to another becomes com- 
parable with the frequency , which characterizes the retardation of the 
electron-phonon interaction. In conventional superconductors the value of 
the coherence length is about 10~^cm and is much lager than a^) ~ lOOA. 
But in high-Tc metalloxide compounds we have the situation with compa- 
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rable with an. Thus, in high- Tc compounds the critical current of the contact 
with dimensions ~ ~ will be sensitive to the effects of strong coupling. 



2 Model and basic equations 

We consider the model of a contact in the form of a filament (narrow channel) 
that joins two superconducting half-spaces (massive banks) (Fig.l). The 
length L and the diameter d of the channel are assumed to be large as 
compared with the Fermi wavelength A^, so we can apply the quasi-classical 
approximation. In the ballistic case, we proceed from the quasi-classical 
Eilenberger equation for the energy-integrated Green's function [^: 

VF.|-G+K3 + A,G'] = 0, (4) 
or 

where 

is the matrix Green function, which depends on the Matsubara frequency uj, 
the electron velocity on the Fermi surface vf and the spatial variable r; 

A(r) = 

is the superconducting pair potential; fj (z = 1, 2, 3) are Pauli matrices. 
Equation for matrix Green's function (^ is supplemented by the normaliza- 





tion condition 11 



= 1. (5) 

The off-diagonal potential A(r) must be determined from the self-consistency 
equation: 

A(r) = A27rT5: </> (6) 

in which < ... > stands for averaging over directions of vp on the Fermi 
surface and A is the electron-phonon coupling constant. In the BCS model 
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the summation over uj contains the cutoff on the frequency cud, which is of 
the order of the Debye frequency . 

The equations and are supplemented by the values of Green's 
functions and A in the bulk superconductors 5*1 and 5*2 far from the channel 
ends : 

Gi,2 = "^^'"1^^''' , Ai,2 = Ao(cos(v9/2)fi ± sin(^/2)f2), (7) 

thus the phase is the total phase difference on the contact. Also we have to 
determine the boundary conditions concerning the reflection of the electrons 
from the surface of the superconductors r^. For simplicity we will assume 
that at r5 electrons undergo the specular reflection. Then for quasiclassical 
Green's function we have the boundary condition (Ref.f^): 

G(vF,r5) = G(vF',r5) (8) 

in which and vp' are the velocities of the incident and specular reflected 
electron. These velocities are related by the conditions, which conserve the 
component of Vi? parallel to the reflecting surface xs and changes the sign of 
the normal component. 

The solutions of the equations (^) and allow us to calculate the current 
density j: 

j(r) = -4z7reiV(0)T < Y^g^ > . (9) 

uj>0 

In the case of the microconstriction shown in Fig.l, under the conditions 
d <^ ^0 and L ^ d {d is the contact diameter), inside the filament we can 
solve the one- dimensional Eilenberger equations with A = A(2;) . The banks 
of the bridge are equivalent here to certain boundary conditions for Green's 
function G{vz,z) at points z = ±-|. Following the procedure, which was 
described in Refs. one can find the Green's functions at the end points 
{z = ±-|) from the general solutions of Eq.(^) in superconducting half-spaces 
Si and 5*2 with conditions (5). They are given by 

G{z = ) = Gi,2 + A,2[Aor3 - (wcos(v?/2) +ir/fisin((^/2))fiT 
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where f2 = ytu^ + Aq^, f] = sign{vz)- The arbitrary constants Ai^2 must be 
determined by matching of these boundary conditions with the solution for 
G{vzz) inside the channel. 

Taking the off-diagonal components of Eq.(|[), we have the following first- 
order differential equations for the anomalous Green's functions: 



^Z I — — V-/i(W5 ^^^^ 



The normal Green's function g^, as follows from condition (|^), is expressed 
in terms of f^j and /J" : 



g. = ^i- (12) 

From equations (j^, (^,(|TT]) and |T^one obtains the symmetry relations 

f:^ivz,z) = {U-Vz,z)r,A*{z) = A{-z) (13) 
and the current conservation inside the channel dj/dz = . 

3 Josephson current and order parameter dis- 
tribution in superconducting microchannel 

In present paper we consider the case of temperatures T close to the critical 
temperature Tc. Near the phase transition curve the order parameter Ao(T) 
in the banks is small. In order to find the Josephson current in the lowest 
order on Aq we linearize the equations (|11] ) on A and obtain ~ Aq(T), 
5f(^ ^ 1 - 1/2/^/+ ~ 1 - O(A^), j ~ A^. The equation for near takes 
the form 

v,^ + 2uU = 2A{z), (14) 
with linearized boundary conditions (|TU| ) 

f^(v.>0,z = -L/2) = ^e-^^, 

Uv.<0,z = +L/2) = ^e^^^. ^'^^ 
Its solution for arbitrary function A{z) is given by 
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z) = ^e-^''¥e-^(^+''^/2)+e-^^ / dz'^^e'^'' . (16) 

to J Vz 

-vL/2 

The Green's function f^{vz, z) is obtained from expression (|14D with the help 
of relations (^ ). 

Substituting function fuj{vz,z) (|l^) in the self-consistency equation (|^), 
we obtain the integral equation for the space-dependent order parameter 
inside the contact 



L/2 

A{z)=A{z)+ J dz'A{z')K{\z-z'\), (17) 

-L/2 

where 



Aiz) = X2nTY. — (e-^^cosH^ + ^^)) , (18) 



K{z) = X2nT (-e-^A . (19) 

The averaging < ... >t,^>o denotes < Fivz = vpcosO) >i,^>o= /q^ c?(cos6')F(cos6'). 

In the case of strongly inhomogeneous microcontact problem the integral 
equation for the order parameter A replaces the differential Ginzburg-Landau 
equation . It contains the needed boundary conditions at the points of contact 
between the filament and the bulk superconductors. Some general properties 
of the solution A{z) of Eq.(^ ) follow from the form of the functions ( 0) 
and (IS). Let us write A{z) in the form 

A(2;) = Ao(T)(cos| + zg(z)sin|) (20) 
and substitute it in the Eq.( |T^). For function q{z) we obtain the equation 

L/2 

q{z)=b{z)+ J dz'q{z')K{\z-z'\), (21) 

-L/2 

with K{z) defined by ( |19]) and the new out-integral function b{z) 
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b{z) = \2nT J2 - (e" - sinh(^)\ . (22) 
In obtaining the Eqs (|2l|) we have used the relation 

X27tT E - = 1' for T T,. (23) 



It follows from (pT]), (p!9| ) and ( p2D that function g(2;) has such properties: 
function q{z) is real, 
g(2;) does not depend on the phase ip, 
in) q{—z) = —q{z),q{0) = 0. 

Thus, the value of the order parameter A at the center of the contact al- 
ways equals to Ao(T) cos ^ . Also, the universal phase dependence of A{z, (p), 
which is determined by (^0|) and i)-iii), leads (see below) to the sinusoidal 
current-phase dependence j = jc sin ip. It is emphasized, that these general 
properties of the ballistic microchannel (within the considered case of "rigid" 
boundary conditions (|10|) and temperatures close to T^) does not depend on 
the contact length L, in particular, on the ratio of L/C,o- 

Now we are going to obtain the Josephson current in the system. To 
calculate the total current I = Sj flowing through the channel at given 
phase difference p>, we use the equation for the current density (^) and the 
obtained above anomalous Green's function /^^ (HI)- The normal Green's 
function ([l2|) in the second order on Ao(T) equals to g^{vz,z) = 1 — 
^fu}{vz, z){f^{—Vz, z))* . It is convenient to calculate the current density at 



the point z = 0. By using the expression for A(z) (pOD , we obtain the general 



formula for the Josephson current I{(p) in terms of function q{z) : 

I{ip) = hsinp,, (24) 



L/2 



v^>0 



dzq{z) 



2uj . 

e ' 



(25) 



Here /q = 7rAQ(T)/(4ei?o7"c) is the critical current at L = 0. It coincides 
with the result of Ref. |§] for the orifice (^ at T near Tc. Expression (|2^ ) 
jointly with equation (|21| ) for g(z) describes the dependence of the critical 



8 



current on the contact length Ic{L). It is vahd for arbitrary value of the 
ratio L/^Q. Note, that in considered here case T Tc, we have the relation 

eo,i^«e(T). 

Let us introduce the dimensionless quantities 

^ = £=![M,^=2n + l,J,= :^. (26) 
Vf rcTc Iq 



In reduced units (Pq ), after taking the average < ... >^^>o, the equations for 
q{x) and Jc take the form 



1/2 

/ 

-1/2 

0][i- 

(2n+l) 



q{x)=b{x)+i J dx'q{x')K{\x -x'\), (27) 



J, = 4 E { <^-Pi-^(2"+i)Hi~^(2»+i)] _ i2Ez[-i{2n + 1)] + 



1/2"=° ' ' (28) 

+Ai J rfa;g(x)[2^£t|gg±iM + 2exEt[-2e{2n + l)x]}, 



where 



b{x) = A E { 2e"PM(2"+i)l^inh[2£(2n+i)x] ^ ^^^^ + 1)(1 - 2x)EzH(2n + 1)(1 " 2x)] + 

n=0 

-£(2n + 1)(1 + 2x)Ei[-£(2n + 1)(1 + 2x)]}, 

(29) 

w 

/s:(a;) = -2A ^ ^i[-2£(2n + l)x]. (30) 

n=0 

Function Ei{x) = •^^p^i^) jg ^j-^g integral exponent. The upper limit N 
in the sums over n is related to the cutoff frequency un in the BCS model, 
~ ujd/Tc. The value of coupling constant A is related to by Eq. (pSf) , or 

N 

in reduced units, 2 A E i2n+i) = 1- the weak coupling limit of A <^ 1, we 

71=0 

have > 1. 

In general case of the arbitrary value of the parameter £ (£ ~ L/^o) the 
Eqs.(^7D are the convenient starting point for the numerical calculation of 
function Jc(^). We consider here two limiting cases of £ ^ 1 and £ ^ 1. 
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For a long microbridge with £ ^ 1 we shall seek a solution of Eq.(|27|) 
in the form q{x) = ax. Substituting this q{x) in Eq. (|27|) , we find that 
a = 2 + 0{l/i). Calculating Jc (^) with q{x) = 2x, we find that the order 
parameter and the critical current are 

A (z) = Ao (^cos I + «^ sin I) , L > ^o, (31) 



/c(L) = — C(3)/o;^,L»^o. (32) 



Expressions (PT|), ( P^D coincide with the solution of GL equations (with ef- 
fective boundary conditions for the order parameter A ) for the clean su- 
perconducting microbridge Thus, our microscopic approach with the 
boundary conditions ( p^Of ) for Green functions (not for A) gives the results of 
the phenomenological theory at L ^ ^q- 

For a short microbridge with £ <^ 1, in zero approximation on i we have 
that q{x) = (A (z) = Aqcos^), Jc = 1. Or, in dimension units, /c(0) = Jq, 
in agreement with formula (|^). The corrections to the zero approximation 
depend on the value of the product iN. For very small i <C Tc/ojd (^-e. 
L <^ an — vf/ujd), the product iN becomes small , although the N ^ 1 . As 
a result, when q{x, £) and Jc{£) are calculated in the region L < an, the cutoff 
in the sums over n must be taken into account. Apparently, when the cutoff 
frequency appears explicitly but not through the value of Tc, the applicability 
of the BCS theory becomes questionable. More rigorous consideration, based 



on the Eliashberg theory of superconductivity |T^, is needed in this case. 
Nevertheless, by using the BCS model with cutoff frequency we suppose 
qualitatively to take into account the retardation effects of electron-phonon 
coupling in our problem. In the domain, defined by the following inequalities 
: < 1, > 1, £ < 1 , the functions b{x) (^ and K{x) (§Up have the 
asymptotes: 

1 1 -\- 2r 

b{x) = 4:\iN{x\n{iN) + x{C + ln2) + — ) + 2xln(l -4x2)]}, (33) 

4 1 — 2x 

A:(|x|) = -2AA^[ln(2£A^ |x|) - 1]. (34) 

Where C ~ 0.577 is the Euler constant. As it follows from Eqs.(^), (|34D in 
this case the integral term in the equation ( pTj) is small, and calculating the 
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critical current in the first approximation on the small parameter IN we can 
put q{x) = b{x). As a result we have 



A{z) = Ao(T)(cos I + ib{z/L) sin |), L < an, 



(35) 



with b{x) defined by expression ), 



Ic{L)=Io{l 



8 T,L 
vrA vp 



),L < an- 



(36) 



In the region {£ ^ 1 and £A^^1} the integral term in the equation (^) 
is numerically small as compared with the out-integral term b{x). By using 
in equation (^) the q{x) = b{x) as the rough approximation, we calculate 
the function Jc{i) shown in Fig. 2. 

For the case £ -C 1, and £N ^ 1, we can put N = oo in the equation for 
q{x) and Jc{i)- The corrections to the critical current in this region of the 
length L can be estimated as 



The expressions (p^), (^) and ( |37D describe the dependence of the critical 
current on the contact length in the limiting cases of short and long chan- 
nel. With increasing length L the critical current decreases. For ultrasmall 
L^aD the value of SIc/Iq ~ (l/A)(L/^o) directly depends on the BCS cou- 
pling constant A, and consequently it is sensitive to the effects of the strong 
electron-phonon coupling. 

4 Conclusion 

We have studied the size dependence of the Josephson critical current in bal- 
listic superconducting microbridges. Near the critical temperature Tc, the 
Eilenberger equations have been solved selfconsistently. The closed integral 
equation for the order parameter A (p!7|) and the formula for the critical 
current Ic (^) are derived. Equations (|17|), ( pSf ) are valid for the arbitrary 
microbridge length L in the scale of the coherence length vf/T^. In 

strongly inhomogeneous microcontact geometry they replace the differential 
Ginzburg-Landau equations and can be solved numerically. In the limiting 




(37) 
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cases L ^ and L -C the analytical expressions for A inside the weak 
link and for the Ic (L) are obtained. In Figure 3 the dependence of Ic (L) on 
L is shown schematically. For long microbridge, L ^ C,o, the critical current 
~ 1/L is in the correspondence with the phenomenological consideration. 
The main interest presents the region L'^C.o, where the microscopic theory 
is needed. We have calculated the corrections to the KO theory (0), which 
are connected with the finite value of the contact size. The expression (2) for 
the Josephson current was obtained in Ref.([§]) in zero approximation on the 
contact size. For the L <^ we obtained that 6Ic/Io ~ — ^ In ^ , where Jq 
is the value of the critical current in KO theory. Thus, the corrections to the 
value Jo are small for L -C ^o, but the derivative dIc/dL has the singularity at 
L = 0. This singularity is smeared, if we take into account the finite value of 
the ratio TcjuJi)- For ultra short microchannel, L^ar, ~ vfIuj£, (the dashed 
region in the Fig. 3), the length dependence of the critical current becomes 
SIc/Iq ~ (-^ ^he constant of electron-phonon coupling). In the very 

small microcontacts we have the unique situation, when the disturbance of 
the superconducting order parameter can be localized on the length ao, mak- 
ing essential the effects of retardation of electron-phonon interaction. The 
ballistic flight of electrons through the channel is dynamical process with the 
characteristic frequency ujq ~ vp/L. For L smaller then ao this frequency 
becomes comparable with the Debye frequency uj£,. 

Thus, the critical current Ic for the finite contact's size is smaller then Jq. 
At the same time the normal state resistance of the ballistic microchannel 
does not depend of the length L and remains equal to the Sharvin resistance 
Rq (H). As the result, the value of the product IcRn is not equal nAl/AeTc 
and depends on the contact size. We have considered here the case of the 
quasiclassical situation, L 3> h/pp- In quantum regime, L ~ h/pp, the 
Sharvin resistance Rq in formula (2) is substituted by quantized resistance 



of the contact, as was firstly shown by Beenakker and Houten [14|. It follows 



from our consideration that for such small microcontacts with L'^ap, the 
rigorous calculation of the Josephson current requires to taking into account 
the retardation effects. 
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Figure Captions 

Figure 1. Model oi S — c — S contact as narrow superconducting channel 
in contact with bulk superconductors Si and 82- 

Figure 2. Dependence of the critical current Jc on the contact length 
L for the microbridge (solid line). The coupling constant A = 0.2. For 
the comparison, the dependence Ic (L) for S-N-S contact ( A = inside the 
channel) is shown (dashed hne). 

Figure 3. Dependence of the critical current on the length of the bridge. 
The asymptotic behavior for short and long bridges are shown. The dashed 
region corresponds to the ultrashort microbridge, L^vf/ujd- 



15 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9904378vl 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9904378vl 



This figure "fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9904378vl 



